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ABSTRACT 


The purpose of this thesis is to investigate a proposed 
reciprocal method for measurements within an exponential 
assembly. It is shown that the usual positions of source and 
detector may be reversed, provided that geometric and macro- 
scopic nuclear correspondence of these two elements is ob- 
tained. Reciprocal and conventional measurements were made 
for neutron density under various conditions of loading and 
source-detector orientation. Certain basic pile parameters 
were calculated using elementary pile theory in order to corre- 
late the two methods of measurement. 


All measurements were made on the Massachusetts Institute 
of Technology subcritical graphite pile. Both indium foil and 
BF3 counting methods were employed. Results obtained definite- 
ly confirm the experimental utility and theoretical accuracy 
of reciprocal measurements. Certain worthwhile advantages are 
associated with the proposed technique, and these are examined 
in detail. 


Thesis Supervisor: Theos J. Thompson, Ph.D 


Title: Professor of Nuclear Engineering; 
Director of M.I.T. Reactor 


ii 





: 


ACKNOWLEDGEMENT 


The authors acknowledge with sincere appreciation 
assistance rendered by: | 


Dr. Theos J. Thompson, thesis supervisor, who 
provided the initial impetus and direction along 
with continuing friendly encouragement; and 


Captain Jeffery Lewins, Royal Engineers, who 
originally conceived the ideas underlying this 

thesis and who furnished advice and guidance through- 
out its course. 


iii 





1bdae 
Je 


TABLE OF CONTENTS 


Title Page 

Abstract 

Acknowledgements 

Table of Contents 

List of Illustrations 

Tables of Data 

Nomenclature 

Introduction 

Procedure 

Results 

Conclusions and Recommendations 
Appendices 

ie Tables of Data 

Jil Statistical Evaluation of Results 
Tit. Details of Procedure 


IV. Supplementary Discussion 


V. Multiplication Ratio 


AE Bibliography 


iv 


Page 


Neal 
nis fa 


iv 


vi 


vii 


2 
a. 
By, 
41 
42 
56 
63 
67 
(2 
15 





Figure 


Figure 


Figure 3 


Figure 4 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


Figure 


10 


ii 


LIST OF ILLUSTRATIONS 


Page 
Photograph of Exponential i 
Pile - North Face 
Photograph of Exponential iL 
Pile - East Face 
Drawing of Exponential Pile 8 
Photograph of Source and 14 
Detector Traverse Equipment 
Horizontal Flux Distribution - a4 
Pile Loaded 
Vertical Flux Distribution - oft 
Pile Loaded 
Harmonic Corrections - Pile 29 
Loaded and Unloaded 
Horizontal Flux Distribution - 33 
Pile Unloaded 
Vertical Flux Distribution - 34 
Pile Unloaded 
Cadmium Ratio 70 
Multiplication Ratio 74 





Table 


Table 


Table 


Table 
Table 
Table 
Table 


Table 


Il 


Iil 


ve 
VII 
VIII 


TABLES OF DATA 


Coordinates of Measurement 
Positions 


Summary of Physical Characteristics 
of Exponential Pile 


Summary of Measured Pile 
Parameters 


Indium Foil Data 

BF3 Data - Pile Loaded 

BF 2 Data - Pile Unloaded 

BF 2 Data - Cd Shutters Installed 


Counter Stability Tests 


vals 


43 


45 


46 


=i 
48 
50 
52 
25 





CR 


Keff 


NOMENCLATURE 


Relaxation length of mn th 
harmonic 


1/Da,; Slope of loge(Ntn) vs. z 
Fermi Age 

Diffusion length 

Material buckling 


Neutron flux density corresponding 
to the mn th harmonic 


Counting rate due to thermal 
neutrons 


East-west coordinate for measure- 
ment of pile dimensions 


North-south coordinate of pile 
dimensions 


Vertical coordinate 

Extrapolated length in x direction 
Extrapolated length in y direction 
Extrapolated length in z direction 
Linear extrapolation distance 


Counting rate for bare foils or 
BF3 tube, no external source 


Counting rate, CD covered, no 
source 


Counting rate, bare with external 
source 


Counting rate, Cd covered, with 
source 


NB Ny 
Cd ratio No or No 


Nen (loaded) 
Multiplication ratio Wi (unloaded) 


Effective pile multiplication 


Infinite multiplication factor 


outa 





I, INTRODUCTION 


A. Historical Summary 


The exponential experiment, in its conventional form, is 
perhaps the classic investigation relating to elementary reac- 
tor theory. Methods of measurement employed and experimental 
data obtained from exponential assemblies were of vital Pt yor 
tance during the early period of nuclear reactor development. 
The first exponential experiment was conceived by Enrico 
Fermi and his associates while conducting their well known 
University of Chicago a iadiee 

Continuing use of these experiments for preliminary 
reactor design information and for student instruction attests 
to the fact that exponential methods have a permanent place 
in the field of reactor technology. 

When intended for design data, the exponential assembly 
is constructed with lattice dimensions identical with those 
of a proposed critical reactor. Overall dimensions are con- 
siderably less, being perhaps one-third of those required for 
criticality. Hence an exponential system will not sustain a 
fission chain reaction because of the relatively great neutron 
leakage which exists. It is possible, however, to attain 
steady state flux conditions by means of an external neutron 
source--usually a Pu-Be or Po-Be cylinder or a reactor thermal 


column. The wave equation for a critical system under steady 





state conditions 
Vee + B® = O (1) 


is closely approximated in the exponential, or subcritical, 
assembly. The general solution. to this equation, following 
the development of Glasstone & Edlund (G-1,ch.9) is 


® = 2 2 A nn cos BEX cos *FYsinh Ying (C2) (2) 


If the lattice has a plane source of thermal neutrons, con- 
tributions of harmonic terms to the flux can be neglected. 


The general solution then becomes 
WT rs : 
= Ancos — cos 4¥ sinhy (c-z) (3) 
Thermal flux distribution along any line parallel to the 
z-axis, then, is given by 


@(z) = F sinh y (c-z)=Ce“[1-e 7 | (4) 


If z does not closely approach c, i.e. if the "end correction" 


term in the bracket is about unity, then 


®(z2)= Ce”? (5) 


Sed 


Hence, in the system under consideration, neutron flux in the 
z-direction diminishes exponentially with distance from the 
origin of coordinates. It is this feature which has given 


the name "exponential" to such subcritical assemblies. 





B. Application of Reciprocity Principles 


The theorem of reciprocity, in its various forms, has 
been valuable in several areas of nuclear physics and reactor 
engineering. Notable was its application by Bethe and Beyster 
(B-3, B-4) for measurement of inelastic neutron scattering. A 
number of other instances are described in the literature. 

Dineet investigation of neutron flux density historically 
has been a primary tool of the experimenter. By this means 
the various parameters required for reactor design are ob- 
tained. In an exponential assembly, flux measurements are 
made in a reproducing lattice supplied with neutrons from an 
external source. These methods have been employed for many 
years, and the concept of neutron flux in this configuration 
is well established. 

It is not illogical, however, that the original concept 
equally well might have been in terms of the adjoint flux 
or importance function. This function is so named because 
the total progeny which result from a single neutron intro- 
duced at a given point in a reactor is a measure of the 
"importance" of this neutron in sustaining the chain reaction. 
Simply speaking, a neutron introduced at the center of a 
reactor, along with its daughter neutrons, is less likely to 
leak from the multiplying region than is the case for a 
corresponding neutron introduced near an edge of the reactor. 
In the case of a thermal, one-group reactor--and only in this 
case--the importance concept predicts an adjoint flux distri- 
bution which corresponds exactly to the flux distribution of 

Be 





thermal neutrons. Consequently, in a static, one-group assembly, 
neutron flux and adjoint flux become identical functions of 
position. 

If consideration is restricted to the one-group reactor, 
the adjoint concept becomes equally acceptable with the con- 
ventional idea of neutron flux, after a brief period of 
orientation. Apart from its wider application in the field 
of multi-group theory, preference is accorded to the classical 
concept because of general familiarity and a well established 
mode of thought. However, if the experimenter is thinking in 
terms of the importance function, certain modifications to 
standard experimental procedures become plausible. 

This thesis is concerned with one of these possible 
modifications which has been proposed by Lewins (L-3). In 
accordance with his suggested application of reciprocity, 
the classical exponential experiment is peered A neutron 
Source is placed within the loaded lattice, as would be done 
intuitively if the experimenter were thinking in terms of 
adjoint flux. A BF3 neutron detector is positioned within 
the preemie pedestal where neutron density becomes a direct 
measure of importance function at the reciprocal source loca- 
tion. In a broad sense it is the purpose of this thesis to 
investigate all facets of such a reciprocal arrangement, with 
emphasis upon accuracy, precision, and experimental utility. 
Theory indicates that reciprocal methods will give results 
equally valid with those obtained by conventional means. It 
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suggests that certain advantages will be associated with 


reciprocal techniques. 


C. The Conventional Exponential Experiment 


In the usual subcritical experiment, detectors are 
placed within a uniformly loaded lattice in positions suffi- 
cient to define the vertical and horizontal distribution of 
neutron density. These detectors typically consist of indium 
foils or boron fluoride gas-filled counters. Since the thermal 
flux distribution is often of primary interest, additional 
measurements may be taken under cadmium covers to separate 
the thermal and epithermal components. It i8 generally ac- 
knowledged that foil measurements provide greater experimental 
accuracy than do BF3 data. This is due largely to the fact 
that foils cause a much smaller local flux perturbation. Under 
certain circumstances, however, experimental results obtained 
from BF3 measurements may be equally acceptable. 

Primary neutrons are supplied alternatively by a reactor 
thermal column or self-contained sources producing neutrons 
by an (@-n) or (yen) reaction. The former method possesses 
a distinct advantage in that the neutrons, as supplied, are 
completely thermalized and are essentially in the foritvof a 
plane source. Use of self-contained sources is generally 
simpler, but harmonic and slowing down effects must be care- 
fully. considered. 

Within limits imposed by pile characteristics and equip- 


ment available, a number of experimental measurements may be 








accomplished on the exponential pile. A partial list includes: 
material buckling (C-2, W-3), diffusion length (R-1, R-2), 

Fermi Age (H-3), effective size (C-3), thermal utilization (R-2), 
effective multiplication (C-2), neutron temperature (H-4), 

albedo (G-1), and reflector savings. These quantities are 
essential in determination of criticality conditions for a 
reactor but are See remeay difficult to determine with suffi- 
clent accuracy from theoretical calculations alone. Further- 
more, the basic exponential experiments enable a student to 
acquire a feeling for reactor physics which is vital to the 


nuclear engineer. 


D. The M.I.T. Exponential Pile 
Photographs of the M.1.T. exponential assembly are 
shown in figures 1 and 2. A detailed drawing appears in 


figure 3. Overall dimensions are as follows: 


x-dimension (east-west) 90.75" (230.5 cm 
y-dimension (north-south) 91.00" (231.1 cm 
z-dimension (vertical) 116.00" (294.7 cm) overall 


25.25" 64.2 cm) pedestal 
90.75" (230.5 cm) lattice 


For uranium loading, 144 fuel channels are provided in 
the lattice region, each channel 1.25" in diameter and extending 
full width in the y-direction. Natural uranium fuel slugs are 
standard AEC issue for educational use--1.087"0D and 8.375" 
long. Nine slugs are loaded in each channel, providing an 
essentially uniform loading on a 7.25" (18.415 em) lattice 
pitch. The assembly is fitted with sliding graphite stringers 


on centerline and extending in both the x- and y-directions. 


Sy 








EAST FACE OF PILE, SHOWING SOURCE POSITIONING 


ROD IN DETECTOR CHANNEL. 


FIGURE 2. 


SOURCE 


NORTH FACE OF EXPONENTIAL ASSEMBLY. 


FIGURE |. 


BLOCKS CAN BE SEEN IN PIT, LOWER RIGHT. 
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Measurements within the lattice may be made with foils or 
BF3 counters. The details of pile dimensions and the coor- 
dinates of measuring positions are listed in Table I of 
Appendix I. A summary of all pertinent physical character- 
istics of the pile is included in Table ITI. 

There are slots at the lattice-pedestal interface for 
insertion of cadmium shutters under those conditions when 
it is desired to minimize passage of thermal neutrons between 
these regions. Shutters for this purpose were fabricated as 
a pyamcadee of this thesis. 

The M.I.T. exponential assembly is physically far 
removed from the reactor and is not designed for utilization 
of a thermal column source. A single level of source channels 
is provided in the pedestal, giving two degrees of freedom 
in source position. A diamond source array, which will mini- 
mize harmonic effects, is possible in the conventional pile 
configuration, though such an arrangement was not used for 


the present experiment. 


E. Purposes of the Thesis 
It was the primary purpose of this thesis to test the 


validity of the theorem of reciprocity as applied to exponen- 
tial measurements. To do this it was first necessary to de- 
vise experimental methods which would permit these measure- 
ments to be made in the reciprocal sense. Following this, 
certain classic pile parameters were to be measured in the 


M.I.T. assembly by both conventional and reciprocal methods, 


Se 





permitting a direct comparison. This comparison was to in- 
clude such factors as statistical accuracy, validity and pre- 
cision of results, and ease of experimentation. 

Particular attention was given to the physical configura- 
tion of source and detector in order that the nuclear enaeae eee 
istics of these components might correspond as nearly as 
possible. Having obtained this correspondence, one could 
logically expect similar perturbations of the flux distribu- 
tion by both source and detector. | 

The counter assembly was examined for anisotropy in 
order to validate measurements taken near the pile boundary 
during horizontal traverses. To accomplish this, a dtirection- 
al stream of fast neutrons was obtained and the counter placed 
in this stream. By taking a series of measurements at 
various angles between beam direction and counter axis, it 
was determined that the degree of anisotropy inherent in the 
BF3 tube was not sufficient to affect measurements within the 
pile. 

Since the reciprocal technique involves much more handling 
of ‘the neutron source than is necessary for conventional exper- 
iments, it was important to devise a safe yet practical means 
for positioning the source. In this connection, radiation 
safety comprised an integral phase of the investigation. 

There are certain results to be anticipated from reci- 
procal measurements beyond those relating to experimental 
accuracy and ease of dnb Bet meeia tion: In addition to verifi- 
cation of the one-group reciprocity within a subcritical pile, 
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a rational, easily visualized substance might be given to 
the concept of adjoint flux. 
With such purposes in mind the following investigation 


was made. 
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II, PROCEDURE 


A. Indium Foil Measurements 

The intended purpose of this phase of experimentation 
was to establish "control" values for pile parameters such 
as buckling and diffusion length, against which BF3 results 
might be compared, Foil data were notably unsuccessful in 
serving this purpose, due largely to low activations and a 
consequent lack of statistical reliability. A complete 
discussion of equipment, procedure, and results is included 


in Appendix C. 


B. Problems Unique to the Reciprocal Me thod 


Two basic difficulties arise in reciprocal measurement 
which are not in common with the usual type of exponential 
experiment. The first of these is associated with the effects 
of source and BF. detector upon pile flux; the second is 
concerned with a possible radiation safety hazard incurred 
by frequent handling of the neutron source. 

In developing experimental procedures for this investi- 
gation, several means for attaining source and detector corres- 
pondence were considered. Generally speaking, these fell 
into one of two categories. On the one hand, differences 
in neutron absorption might be minimized by fitting both 
source and detector within suitable containers of large macro- 
scopic absorption. This method must accept very large pertur- 
bations of the flux and consequently must sacrifice a considera- 
ble degree of accuracy in defining its distribution. A second 


ve 





approach to this problem consists of matching certain un- 
avoidable differences between source and detector as closely 
as possible, accepting whatever variance remains. This latter 
method was ultimately employed. 

Two Pu-Be Sources were used, with a combined overall 
length closely approximating that of the active portion of 
the BF3 tube. A source container was fabricated in the shop 
from aluminum tubing. Since the container is thin and the 
neutron cross section for aluminum small, it is plausible to 
neglect any effect of this component upon the neutron flux 
distribution. To this container was attached a base assembly, 
complete with coaxial cable, from an unusable BF 2 tube of the 
Same type as that employed for counting. For epithermal 
measurements, suitable cadmium jackets of identical weight 
were made for both source and detector. Figure 4 is a photo- 
graph of these assemblies as they finally evolved. 

In this manner good correspondence was obtained, assuming 
the total macroscopic cross section of source cylinders and 
gas-filled portion of the tube to be eeauiee |. This assump- 
tion, though by no means precise, is not a bad one. Calcula- 
tions were made to determine the total macroscopic cross 
sections of the active portion of the BF 3 tube and of the 
neutron sources. The computations showed that the primary 
neutron absorbers, namely BLO and Pue39, had total neutron 
absorption cross sections of the order of .6 cm and 8 cm@ 
respectively. it is felt that the desirably small flux per- 
turbation which arises is adequate to justify a slight lack of 
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correspondence between source and detector. 

The second basic difficulty concerns personnel safety 
during source handling. Since the source must be moved once 
for each separate positioning of the detector in a conventional 
experiment, it was necessary to devise an effective yet 
simple method of precise control which would involve a mini- 
mum of exposure to the source. To accomplish this, a sec- 
tion of 3/8" aluminum rod approximately six feet long was 
connected to the source container. A corresponding length 
was attached to the detector to insure similar effects upon 
flux distribution. These handling rods afforded a means of 
positioning source and detector from the pile surface, pro- 
viding excellent accuracy of position. It is estimated that 
all measurement locations are correct to 1/32" or less. Hence 
the possible effect of positioning error is not considered 
in evaluating the results of this investigation. 

Film badges and pocket dosimeters were worn by both 
thesis participants throughout all experimental phases, thus 


providing a final evaluation of radiation safety aspects. 
c.. BF 2 Measurements -=- Pile Loaded 


This phase of experimental work consisted of eight com- 
plete sets of vertical and horizontal traverses within the 
reproducing lattice. Four series of measurements were made 
which corresponded exactly to the classic experiment for 


material buckling. These consisted of BF 3 traverses with the 
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detector alternately bare and shielded by cadmium for both 
source and no-source conditions. Measurements were made at 
each of seventeen different positions within the assembly, 
eleven of these defining the vertical flux and six being 
associated with the horizontal distribution. 

Following this, another series of four measurements was 
made in the reciprocal configuration. Two sets consisted of 
bare and cadmium covered source traverses throughout the 
seventeen lattice positions. The detector remained in a 
fixed pedestal location, being alternately bare and sheathed 
with cadmium. Measurement of reciprocal “background” was 
accomplished simply by two counting periods -- bare and 
cadmium -= with source removed from the lattice and detector 
stationary within the pedestal. 

A general outline of the conditions under which BF3 
counts were made in the loaded pile, with reference to 
appropriate data tables in Appendix I, is given here; sta- 
tions referred to are shown clearly in Figure 3, and their 
coordinates are given in Table I, Appendix I. 

Run 1 Source in pedestal: count with detector at 

each of stations 1 - 11 anda - with and 
without cadmium covers in place. (34 counts) 

Run 2 No source; same procedure as above, to deter- 
mine thermal neutron background from spontaneous 
fission. (34 counts) 

Run 3 Detector in pedestal: count with source at 


each of stations 1 - 11 and a - g, with and 
without cadmium covers in place. (34 counts) 
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Run 4 Detector in pedestal: count with no source 
in lattice, with and without cadmium covers in 
place, to determine the spontaneous fission 
thermal neutron background. (2 counts) 
(Data from runs 1 - 4 appear in Table V, Appendix I). 
As a general policy, duration of counting was adjusted 
for each position in order that final statistical error 
Should not exceed 1%. Data were corrected for resolving 
time of the BF3 tube whenever these corrections approached 
statistical significance (E-1). 
Conventional and reciprocal data were adjusted for 
end leakage and diffusion harmonics by iterative methods. 
After considerable thought it was determined that these 
corrections are identical for corresponding conventional 
or reciprocal measurements. This is equivalent to stating 
that the corrections are functions of position and relaxa- 
tion length alone, applying equally to the one-group distri- 
butions of flux and adjoint flux. The correction for 
diffusion harmonics was applied to a two-region system, employ- 
ing different harmonic attenuation within pedestal and lattice 
portions of the pile. The basic procedure employed is de- 
tailed by Glasstone and Edlund (G-1,p.126). <A modification 
was introduced by separate consideration of harmonic attenu- 
ation with pedestal and lattice portions of the pile. Since 
these two regions are of totally different character when 
the lattice is loaded with uranium, it is evident that the 
relaxation length for a given harmonic must change abruptly 


at the interface. 
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Assuming the pedestal to consist of solid graphite, 
K*® from the unloaded data was adjusted for density varia- 


tion as later shown by equation (13) and Ymn calculated from 
a z mM \pe n® \e 
van = & iF [( a ) + ( b ) ] (6) 


Thence a harmonic correction factor may be obtained from 


the relationship 


Ge = | + ee e~ %ma X72 | a) 


n 
and the fundamental mode Sireaied by division of the 
measured density by this factor. By substitution of the 
appropriate value for z in equation (7), the magnitude of 
Cy at the interface is obtained. 


Within the reproducing region, Ym, is derived from 
ot. mz \* nv7,*7 2 
yc = [(@F) + (¥) ]- B, (8) 
and a similar procedure employed to calculate Cy at each of 
the vertical measurement positions. 

All of these calculations are iterative in that the 
higher harmonics are first neglected in calculating y from 
(6) or (8). Preliminary values for Cy at each measurement 
position are thus obtained and applied to the data. This in 
turn yields a new Y, and the procedure is repeated until 


successive approximations cause no change in the slope. It 


is necessary to conduct these calculations concurrently with 
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the iteration for end leakage. It was found that accuracy 

within the limits of counting statistics was generally ob- 

tained by (1) fitting the slope by eye (2) applying correc- 
tions (3) fitting altered slope by eye (4) applying correc- 
tions, and finally (5) fitting the slope by the method of 

weighted least-squares. 

The same equipment was used for all BF3 measurements 
reported in this thesis. High voltage was applied to the 
detector tube approximately 48 hours before any counts were 
made. Counter-scaler settings were not changed during the 
entire experiment. Data were not corrected for deviation 
from standard temperature and eit s dine: since such correc- 
tions are variable and extremely small. An adequate check 
upon counting stability was provided by periodic measurements 
at several standard positions. No statistically valid 
trend could be observed; hence no stability compensation 
was applied. The data are reported in Table VIII of 
Appendix I. 

The neutron sources used are 2.54 x 3.65 cm Pu-Be 
cylinders clad in tantalum and Stainless steel. Each of 
the two sources contains one curie of Pu©29 and emits 
approximately 1.6 x 10 fast neutrons per second. Average 
energy of the source spectrum is 4 Mev (S-3). The chamber 
of the Nancy-Wood BF. detector used is 10 cm long, 1.59 cm 
in diameter, has an active length of 8.25 cm, and is filled 
to a pressure of 40 mm of Hg with BF 3 enriched in boron-10. 
This counter was connected directly to a binary scaling circuit. 
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DD BF Measurements -= Pile Unloaded 


After all uranium had been removed from the pile, a 
further series of measurements was made to determine neutron 
diffusion length in both the conventional and reciprocal 
arrangements. First the usual sigma pile experiment was 
conducted, with source in the pedestal and detector tra- 
verses made through the same seventeen stations in the un- 
loaded lattice. This again included both bare and cadmium 
covered measurements. Pile background was very low and 
essentially independent of position; thus it was found 
that single background counts with and without cadmium 
cover, were sufficient to define cheneal meltrou background 
throughout the unloaded pile. | 

Reciprocal measurement was accomplished by reversing 
source and detector and repeating the conventional traverses. 
Background was particularly low in this configuration, since 
the detector remains within a well shielded pedestal. 

Again an outline of conditions under which BF 3 counts 
were made in the unloaded pile is given; Figure 3 and Table I 
of Appendix I specify stations. (Discontinuity of run 
numbers is due to the fact that runs were numbered chrono- 
logically). 

Run 12 Source in pedestal: count with detector at 
stations 1 - 11 and a = g, with and without 
cadmium covers in place. (34 counts) 

Run 13 No source in pile: count with detector at 
stations 1 - 11 and a = g, with and without 
cadmium covers in place, to determine thermal 
neutron background. (34 counts - found sta- 


tistically constant) 
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Run 14 Detector in pedestal: count with source in 
Stations 1 - 11 and a - g, with and without 
cadmium covers in place. (34 counts) 

Run 15 Detector in pedestal: background count with 
no source in pile, with and without cadmium 
cover in place. (2 counts) 

All comments of section II. C. regarding statistics, 

corrections, and equipment apply equally to these unloaded 


measurements. 
E. BF 3 Measurements - Shutters Installed 


The BF data mentioned during discussion of loaded and 
unloaded measurements were designed to provide experimental 
values of pile buckling and diffusion length. Moreover, the 
methods employed furnished a direct evaluation of reciprocal 
methods. 

It was felt that additional measurements might provide 
a better insight into the importance function. Specifically, 
flux determinations were made with both source and detector 
within the loaded lattice. Cadmium shutters were fabricated 
from sheet stock and inserted into slots provided at the 
lattice-pedestal interface. Reflection of thermal neutrons 
from the pedestal was thus minimized, and lattice boundary 
conditions were made more nearly symmetric. A total of 
twelve complete traverses were conducted to measure thermal 
and epithermal neutron density. These consisted of alter- 
nately positioning bare or cadmium covered source and 
detector about the central and two extreme vertical measure- 
ment positions. 
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BF, counts made with shutters installed are here out- 


lined, specifying conditions of measurement. 


Run 


Run 


Run 


Run 


Run 


Run 


Run 


-) 


10 


11 


Cadmium shutters installed, source in station 
6(f): count with detector at stations l - 5 
and 7 - 11, with and without cadmium covers 
in place, (20 counts) 


Cadmium shutters installed, detector in 
station 6(f): count with source in stations 
1 = 5 and 7 - 11, with and without cadmium 
covers in place, (20 counts) 


Cadmium shutters installed, no source in 
assembly: count with detector in each of 
stations 1 - 11 and a - g, with and without 
cadmium covers in place. (34 counts). This 
is similar to Run 2, but with cadmium shutters 
installed. 


Cadmium shutters installed, source in station 
ls count with detector at stations 2 - 11 and 
a- g, with and without cadmium covers in 
place. (32 counts) 


Cadmium shutters installed, detector in station 
l: count with source in stations 2 - 11 and 
a- g, with and without cadmium covers in 
place, (32 counts) 


Cadmium shutters installed, source in station 
lls count with detector in stations 1 - 10 and 
a- 2, with and without cadmium covers in 
place. (32 counts) 


Cadmium shutters installed, detector in station 
ll: count with source in stations 1 - 10 and 
a- 2g, with and without cadmium covers in 
place. (32 counts) 


Data from runs 5 - 11 appear in Table VII, Appendix I 


Previous mention of statistical accuracy and resolving 


time is applicable to these data. Further corrections were 


not applied. 
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III, RESULTS 


A. Foil Measurements 

Data obtained through indium foil measurements did 
not contribute to the information required of this experi- 
ment. Statistical errors prohibited accurate calculation 
of material buckling of diffusion length. This factor com- 
bined with insufficient data in the slowing down region to 
prevent experimental determination of Fermi Age. Consequent- 
ly the data are reported in Table IV, but no further cal- 


culations are made. 
B. BF 3 Measurements -- Pile Loaded 


The primary purpose of these measurements was to pro- 
vide experimental values for material buckling obtained by 
conventional and reciprocal means. Correlation of results 
for these two methods was excellent, strongly indicating 
the validity of reciprocal data. Certain very definite 
advantages are associated with the novel procedures pro- 
posed in this thesis, particularly with regard to statisti- 
cal accuracy and precision. All data for these measure- 


ments are included in Table V. 


1. Horizontal Distribution 

Figure 5 shows the horizontal flux distribution at 
mid-height under both conventional and reciprocal condi- 
tions. Plotted points indicate experimental results 


and standard error in the data. These are referenced 
23. 
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FIGURE 5 


HORIZONTAL DISTRIBUTION OF 
THERMAL NEUTRON DENSITY AT 
LEVEL 6 (1565 CM ABOVE SOURCE 
LEVEL) IN THE LOADED ASSEMBLY. 
LINES SHOW BEST COSINE FIT. 
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to a cosine distribution, which was normalized by the 
method of least-squares. A reasonably good cosine fit 
was obtained in both cases. Mathematical fitting of 
these curves was undertaken, since they were used quan= 
titatively to deduce the horizontal linear extrapo- 
lation distance. Both conventional and reciprocal 
curves converge to an extrapolated boundary 0.75" 

(1.9 cm) beyond the physical dimension, which is char- 
acteristic of other experimental values (K-1). 

This value was also used in calculations involving 
the vertical extrapolation distance. More properly, 
the extrapolated height should be obtained by trial 
and error, minimizing squared deviations of the experi- 
mental points for different values of this dimension. 
This is prohibitively tedious if manual methods of 
calculation must be used, since the process becomes 
doubly iterative with the fitting for end and harmonic 
corrections. | 

It should be noted that any reasonable curve 
drawn through the conventional data has a pronounced 
dip or flux depression between stations "g" and "e", 
centered approximately on the pile centerline. Pile 
loading was uniform in this region and there are no 
geometric irregularities in the pile which might con- 
tribute to this effect. The phenomenon was observed 
during loaded and unloaded traverses made in the con- 


ventional manner and has not been satisfactorily explained. 
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2. Vertical Distribution 

Results of conventional and reciprocal traverses 
made in the vertical direction are shown in Figure 6. 
All data points were corrected for diffusion harmonics and 
end leakage. Standard error is indicated by brackets on 
the individual points. Fitting was by the method of 
least-squares for stations 3 through 11, points being 
weighted inversely with their variance. The slopes 
of these lines, and hence the pile bucklings, are iden- 
tical within statistical error. 

Certain anomalies are apparent in the results, and 
these are discussed more thoroughly in Appendix IV. A 
qualitative physical explanation for the large differ- 
ence in counting rate has been developed. Furthermore, 
as may be seen from Figure 6, conventional data for 
stations where source and detector are in closest 
proximity exhibit an appreciable effect due to slowing 
down density and change in energy distribution. This 
effect is absent in the reciprocal configuration, where 
inclusion of Stations 1 and 2 in a weighted least-squares 
analysis does not alter the results. Appendix II shows 
that these two reciprocal measurements are statistically 
acceptable. 

Material buckling was experimentally obtained by 
each method. These values are; 

Conventional Reciprocal 
Bh = 86412 x 107-8 em-2 Ba = 80#11 x 10-68 em-2 
26n 
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FIGURE 6 


DETERMINATION OF T.. IN NATURAL URANIUM — GRAPHITE 


LATTICE 8Y CONVENTIONAL EXPONENTIAL EXPERIMENT 


AND BY THE RECIPROCAL METHOD, WITH SOURCE AND 
DETECTOR INTERCHANGED. 
STATIONS 3-11, WITH END CORRECTIONS AND HARMONIC 


CORRECTIONS (THROUGH 5,5) WERE USED IN WEIGHTED 
LEAST - SQUARES ANALYSES TO FIND Tu AND BUCKLING. 
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A superficial inspection would indicate that the 
standard error reported for these results is excessive- 
ly high, considering the 1% accuracy in counting statistics 
which was maintained throughout the experimental work. 
The error reported is appropriate,resulting in part 
from the propagation of counting error but largely 

from a small difference in the subtraction of large 
numbers. This may be readily visualized by examination 
of equation (8), from which the results are obtained. 
Valid comparison with other experimenters is difficult 
because of differences in lattice arrangement. For 
lattices of similar design, material bucklings of 

about 55 x 10-6 to 88 x 10-6 em7= have been obtained 
(D-1, G-5). 


3. Corrections to Data 
Referring to the development on pages 17 and 18, 
experimental points were corrected to the fundamental 


mode by the relation 


_ _Ntp(measured) 
N,,( fundamental ) = C,, Ge (9) 


where end leakage is compensated by the expression 
-27 (¢c- 
Ce.= |- e *? e-2) (20) 
and the effect of diffusion harmonics is corrected by 


Cys 1+ [Sb emi | (7) 
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See FIGURE 7. 
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where 

Yon = [ (82) + (ME) J- Be (8) 
In these calculations, Y, is the reciprocal relaxation 
distance of the fundamental harmonic; Br is the material 
buckling. The above expressions are standard for such 
corrections. Harmonic corrections were made on a two 
region basis through the 5th harmonic. This consisted 


of calculating an attenuation of each harmonic within 


the pedestal region from (6) employing 
2 2 2 2 
Ym = kK + [(GT) + (SE) J (6) 


Having the relative harmonic strengths at the lattice- 
pedestal interface, it was a relatively easy matter 

to attenuate these harmonics exponentially according 
to equations (6) and (7). This was done graphically 
on the semi-logarithmic plot of Figure 7, saving much 
labor of calculation and sacrificing no meaningful 
accuracy. It should be noted that end and harmonic 
corrections are by no means negligible at any position 


in a small assembly employing a point source of neutrons. 


4, Evaluation of Technique 


The statistical benefits of reciprocal measurements 
have already been mentioned and were graphically shown 
in Figures 5 and 6. For any desired degree of accuracy, 


actual counting times in a given position are reduced 
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over 25% by virtue of the greater counting rates ob- 
tained in the reciprocal configuration. The precision 
of reciprocal measurements is far superior, as discussed 
in Appendix II. 

Handling of the source can be accomplished safely 
and effectively without elaborate equipment, though 
the procedure requires respect and appreciation for 
possible danger involved. During the four week period 
of thesis experimentation, ionizing radiation received 
by the participants was never significantly above back- 
ground. The average accumulated neutron dosage was 20 
millirem, as measured by tracks counted on a fast neutron 
film. This dosage, again, was not significant. 

An absence of detector handling, on the other hand, 
presents decided advantages. Permitting the BF tube to 
remain in a single position throughout an experiment 
promises better electronic stability, less likelihood of 
damage to the tube, and a probability of more consistent 
results. 

Calculation procedures are considerably apbrewiaed 
for reciprocal data, since pedestal background is 7 
measured once and remains constant for all source posi- 
tions. Since fewer measurements are required and higher 
counting rates are obtained in the reciprocal arrange- 
ment, the overall time required for experimental work 


is approximately 40% that of a conventional experiment. 
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C.. BF3 Measurements = Pile Unloaded 





After all fuel was removed from the exponential assembly, 
measurements were made to permit determination of neutron 
diffusion length in the graphite used in this assembly. These 
data were obtained by conventional and reciprocal awe 
ments, and correlation of results was again excellent. Although 
reciprocity was thus shown to be valid in an unloaded assembly, 
certain of the advantages applying to the reciprocal method 
in an exponential experiment were of less consequence in this 


phase of the investigation. 


1. Horizontal Distribution 

Horizontal measurements of neutron density by 
reciprocal and conventional methods are shown in Figure 
8. These data permit an acceptable fit with an approxi- 
mately normalized cosine distribution. Linear extra- 
polation distance for horizontal flux was determined 
to be 0.75" (1.9 cm) -- identical to the result obtained 
in a loaded lattice. This value was again used for verti- 
cal extrapolation, for reasons previously discussed. The 
conventional data indicate a central flux depression 
Similar to that observed during measurements performed 


in the loaded pile. 


2. Vertical Distribution. 
Figure 9 includes the results of conventional and 


reciprocal traverses along the z-axis. Measurements 
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DISTRIBUTION OF THERMAL 
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were correctéd for end effect in the usual manner. Two 
region harmonic corrections were again applied as pre- 
viously outlined. The need for a two region analysis 
is nner leas for the unloaded data, since relaxation 
length differs between regions only because of the 
density change due to air channels within the lattice. 


Diffusion length as determined by the two methods 


is: 
Conventional Reciprocal 
L = 46.140.9 cm L #45,540.8 cm 


Corrected to standard (1.6 gm/cm3) graphite, these 


results become: 


Conventional Reciprocal 
L = 48,.7#%1.0 cm L248.1£0.9 cm 


These values are somewhat below a world-consistent 
measurement of 52.0¢1.0 cm (K-1). The discrepancy is 
attributed to impurity of the surplus graphite used 

in construction, and to other impurities which have 
accumulated within the pile. More important to the 
purposes of this thesis is the fact that statistically 
identical results were obtained by the two methods, as 


Shown in Appendix II. 


3. Corrections to Data 
Equations (6), (7), (9) and (10) of Section IV.B.3. 
were used to correct experimental measurements for diffu- 


sion harmonics and end leakage. This process, which 
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again employed graphical methods, requires no further dis- 

cusSion. Figure 7 illustrates harmonic effects. | 
Corrections to compensate for cylindrical voids in 

the moderating medium were investigated. According to 

Wilson, et al (W-3), the ratio of diffusion length 

parallel (Ly) and perpendicular (Ly) to a channel 

axis is given by: 


(Lily = ce[ it 2v(lt ry) eh 


Li’ Left y(2t t7 ae) | 


where WY is the ratio of void volume to moderator 
volume, r is the void hydraulic radius, and is the 
transport mean free path in a homogenized lattice. 


Thence the expression for diffusion length becomes: 


ne = kt = Yu? - (E+E) (12) 
It was determined that this correction is totally 
negligible in the M.I.T. assembly. 

For comparison with published data it was necessary 
to correct the experimental diffusion length to standard 
density, temperature, and pressure. This was done by 


the method presented by Wilson, et al (W-3): 
L*(1.6)= L? [ I+(I/p,- 1/p,) 20-5 1.925-10°* ](P2/1.6)° (13) 


wheres: Psatmospheric pressure in mm of Hg 
T sgraphite temperature in °K 
Po? density of crystalline graphite 
p,zaverage density of the lattice including 
empty channels 
P2zaverage density of solid graphite in the 
lattice 
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Og = weighted average neutron cross section 
of air (1.474 barns) 
O = cross section for pure carbon (.003748 barns) 


It should again be noted that harmonic effects and 
end leakage are of consequence at all positions of measure- 


ment. 


4, Evaluation of Technique 
Figures 8 and 9 indicate a slight advantage in 


counting rate associated with the reciprocal method. 

This has been attributed to geometrical differences 

and is fully discussed in Appendix IV. Benefits of elec- 
tronic stability mentioned in evaluation of exponential 
techniques apply equally to the nmloaded measurements. 

In summary, it may be said that reciprocal methods 

are completely acceptable for the sigma pile experiment 
and that a considerable improvement in precision of 
measurement is gained. The effect of experimental 
technique upon precision is discussed more fully in 


Appendix II. 


S. BF 3 Measurements -=- Pile Loaded, Shutters Installed 





These measurements, originally designed to yield addi- 
tional insight regarding behavior of the importance function, 
consisted of several series of flux traverses under varying 
conditions of source=-detector orientation. In all cases, 
lattice and pedestal were isolated by cadmium shutters in- 


stalled at the interface. Measurements by conventional and 


ae 





reciprocal methods defined the thermal flux employing stations 
1, 6(f), and 11 as reference points. 

The data thus obtained has been corrected to relative 
thermal density and qualitatively examined. The validity 
of reciprocity and the adjoint flux concept was established 
by one-group analyses of experiments conducted without the 
cadmium shutters in place; no reason to modify previous Sea: 
clusions was inferred from the shutter data. 

It is possible that a more thorough analysis might eros 
vide information which would yield better definition of 
boundary conditions, reflector effect, and fast fission. 


Such a study is beyond the scope of this thesis. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 


Wye 


Measurement of pile parameters in an exponential 
assembly containing a single source may be per- 
formed equally well when source and detector are 
interchanged. 

Corrections to reciprocal data for end leakage 
and harmonic effect may be applied in the usual 
manner. The methods for obtaining these correc- 
tions remain unchanged. 

Accuracy of measurement is considerably improved 
by use of reciprocal methods, and somewhat fewer 
measurements are required to define the flux dis- 
tribution. Hence the time required for experi- 
mental work is reduced appreciably. 

Calculations required to reduce reciprocal data 
to thermal density are shortened, diminishing 
the time required for computation of parameters. 
All points of measurement are equally valid in 
the novel experiment. 

The precision of reciprocal measurements is sig- 
nificantly greater than that obtained by conven- 
tional methods. 

No appreciable radiation hazard is associated 
with reciprocal techniques, provided reasonable 
precautions are observed. 
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Information regarding neutron energy distri- 
bution in a loaded lattice may be obtained 
by comparison of conventional and reciprocal 
experiments. This possibility is discussed 
in Appendix IV. 


B. Recommendations 


7. 


Educational and research institutions should 
consider the reciprocal technique + a possible 
improvement in the traditional method of per- 
forming an exponential experiment. A presenta- 
tion of the physical meaning of importance is 
implicit in the novel approach. 

Theoretical calculations and experimental work 
should be performed to verify or disprove the 
validity of measurements made with four detec- 
tors in the pedestal and a single source in the 
lattice. This arrangement should provide can- 
cellation of certain major harmonics. 

Further investigation of the relevance of this 
experiment to definition of neutron spectra is 
recommended. Appendix IV contains some ideas 


in relation to this aspect of the experiment. 
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APPENDICES 
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APPENDIX I 


TABLES OF DATA 


The following tables present all experimental 
data obtained during the experimental phases of this 
thesis. Certain further information is included re- 
lating to physical and nuclear characteristics of the 
M.I.T. exponential assembly. 


All measurements of neutron density have been 


corrected to the extent indicated by notes following 
the individual data. 
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TABLE I 


Coordinates of Foil Measurement Positions 


Position 


11 
10 


OO PNWEUD TO AOD HR OWN OO 


Notes; 


(1) 
(2) 


X y Z 
-4 8125 -1.75 83.25 
~4.8125 76 
-4.8125 68.75 
-4 8125 : 61.50 
-4 8125 54.25 
-4 8125 U7 
-4,8125 U7 
-10.6250 47 
~16.4375 47 
~22,2500 47 
-28 .0625 UT 
- 33.8750 47 
~ 39.6875 47 
-4 8125 39.75 
-4,.8125 32.50 
~-4 8125 25.25 
-4 8125 18 
-4.8125 10.75 
~1.6250 0 ~16.25 


z20O at lattice-pedestal interface 
All coordinates expressed in inches from 


intersection of pile centerline and zs0 
plane. 
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TABLE I (Continued) 


Coordinates of BF Measurement Positions 


Position X y Z 
11 0 =e aio 81.625 
10 0 74. 375 

3 O 67.125 
8 O 59.875 
7 0 52.625 
6 0 45,375 
g -(.250 45.375 
f AS. 375 
e xeno 45.375 
a 14.500 45,375 
Cc 217 150 45.375 
b 2 45.375 
a 36.250 45.375 
5 38.125 
4 0 30.875 
3 0 23.625 
2 0 16. 375 
1 O 9.125 
Ss ~1.625 O ~16,250 


Notes: 
~ (1) zsO at lattice-pedestal interface 


(2) All coordinates expressed in inches from 


intersection of pile centerline and z=20 
plane. 
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*Loading 
loading. 


TABLE II 


SUMMARY OF PHYSICAL CHARACTERISTICS OF 
EXPONENTIAL PILE 


Linear Dimensions 
x 90.75" (east-west) 


y 91.00" (north-south) 
z 90.75" (vertical-lattice only) 


Volumes 
Overall volume, including air channels 749,436 in?> 
Volume of air channels 17,205 in3 
Volume of graphite (Vq) 732,231 in3 
Volume of fuel (Vy) 8,099 in? 
Ya = go.41 
Vu 
Fuel 
Uniform loading 9 slugs per fuel channel 
No. channels 144 
Total slugs* 1289 
Wt. per slug 4.275 1b 
Slug dimensions 1.087" OD 
8.375" length 
atoms/cm3 (Ny) 4.7325 x 1022 
Moderator 
Density 1.662 gm/em3 
atoms/em> (Ng) 8.3370 x 1022 
No/Ny x VG/Vu 159.27 
Ny/Ng x Vu/Ve 006279 


in isolated channels departed from the uniform 
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TABLE III 


SUMMARY OF MEASURED PILE PARAMETERS 


Conventional Reciprocal 
Material Buckling 86412x10-5em"2 80¢11x10~Sem-2 
B, Diffusion Length (unloaded) 46.12.90em 45.5 +.8em 
Diffusion Length (unloaded) 
(corrected to std. conditions) 48.741.0cm 48.14 .9cem 
D. Diffusion Length (loaded)* 14.71 cm 14.54 em 
E, Fermi Age* 354 cme 354 em? 
F, Infinite Multiplication 
Factor (Ko )* 1.049 1.045 


*Calculated value 
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TABLE IV 


Indium Foil Data 


Foil N N N CR N 
an B C b rel 
wat 21. 33 6.71 13.10 3.18 4.97 
10 43,81 13.50 13.27 3.25 20.55 
9 64.61 Vi 17.84 3.69 33.95 
8 104,66 20.17 21.16 5.19 68.90 
7 132.48 28.01 17.23 4.73 91.77 
6 144.59 41.07 20.58 3.52 87 . 36 
g 144,59 41.07 20.58 3.52 87 . 36 
f 147.03 35.81 23.67 4,11 83.79 
e 153.69 28.69 20.69 5 Er 109.76 
d 110. 33 26.40 18.69 .18 70.16 
Cc 96.46 20.76 20.44 4.65 60.63 
b 74,14 15.16 21.81 4.89 4e.89 
a 34.56 Viet 22.06 3,08 7.09 
5 203.28 50.89 19.87 3.99 V3 ieaes 
4 302,40 53.50 23.27 5.65 231.65 
3 434.12 67.56 17.98 6.43 353. 31 
2 604.50 124,50 17. 38 4,86 467.20 
1 937.63 129.13 2 70 7.26 782.42 


These data do not imclude harmonic or end corrections. 
The number of significant figures reported is a result 
of calculation methods employed and is by no means a 
measure of statistical accuracy. 
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TABLE V 


BF'3 Data - Pile Loaded 


OB Nc QO 


6 12 1 4 
ik 18 1 74 

9g at 2 92 

10 38 «62 104 
14 56 2 107 
18 73) i 116 
19 72 3 120 
18 13 Soe dis 
19 69 35 119 
ile’ 65 3112 
16 54 2 102 
12 33 2 78 

8 a4 e2 ug 

23 101 Jaa 
36 139 1 115 
52 191 1 107 
66 285 2 100 
120 399 2 92 


Reciprocal (Detector in pedestal 


Posit 


ak 
10 


Mow Fu oO OO FIR ON OVO 


yo 

756 

1170 
1684 
2 336 
Boles 
3250 
S315 
3120 
2799 
2297 
1723 
10 32 
4579 
6562 
9411 


13,817 
21,277 


5 5 1 3 
9 3 a 
12 8 a 
18 9 il 
22 eS 
a4 20a! 
2S ie); wall 
24 20 1 
23 ie aL 
a4 iS al 
21 11 
19 kale ah 
10 8 1 
48 28 8 
eM 39 e 
97 Do alee 
ei ae 
146 188 4 
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Stim. 2k 
Soe 1 
3 4 |] 
bie gee 
8. 3D. ool 
ec .e2 l 
Sis ee 
PIs apa 
3. 5 ee 
3° Sa 
3° yee 
source in 
oO 0. oO 


Nb Ob Ne Ge 


osition 
Neh Oth Neh 
Corr 
OUT 6 4og 
408 8 552 
639 9 (22 
914 Thal 965 
1308 14 1328 
1842 19 1817 
1890 19 1881 
1842 19 1817 
1883 19 1874 
1687 18 1684 
1422 17 1423 
1034 12 1039 
624 9 6 30 
2529 e4 eho6 
Sys 36 3439 
5170 52 4545 
8404 66 6796 
13,927 120 9881 
osition 
4O7 5 708 
720 9 928 
1129 12 ens 
1642 18 17 34 
2291 22 862325 
3260 ey 3216 
3199 23 3133 
3260 a4 3216 
3069 23 3054 
2751 ey e746 
2253 21 2255 
1679 19 1687 
991 10 1001 
4518 48 =6§4 334 
6490 57 6011 
9319 97 8192 
13,687 117 liee6e 
21,056 146 14,939 





TABLE V 


Notes: 


(1) Resolving time corrections have been applied 
to Ng when appropriate. 


(2) Nn (corrected) includes corrections for 
diffusion harmonics and end leakage 


(3) All data is rounded to the nearest whole 


count. This was done after the best values 
had been finally corrécted. 
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Conventional (Source in pedestal, detector in position) 


Posit 
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BF3 Data - Pile Unloaded 
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TABLE VI 


Notes: 


(1) 


(2) 
(3) 


(4) 
(5) 


Resolving time corrections have been applied 
to Np when appropriate. 


Nin (corrected) includes corrections for 
diffusion harmonics and end leakage 


All data is rounded to the nearest whole 
count. This was done after the best values 
had been finally corrected. 


The standard deviation of No was not sta- 
tistically significant. 


Np, as reported, is corrected for spurious 

background. This was 3.0 cpm in the conven-~ 
tiomal case - 2.0 cpm for reciprocal measure- 
ments. Spurious background for cadmium data 
was 0.1 comventional and zero for reciprocal. 
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TABLE VII 
BF 3 Data - Shutters Installed 


Source Position 1 - Detector Traversed 


Tube 


Np 


Os 


Nc 


Oc 


Np 


ot 


= 
re) 


a 


N 


Posit ve 
11 #4908 10 35 2 52 5 2 1 626 
10 1620 13 58 2 78 6 3 1 1488 
9 2hoo 22 92 3 86 [T 5 2 £2228 
8 3525 34 130 5 ay (G2 dy al 3281 
7 5131 51 192 6 107 7 4 1 4836 
6 7473 43 260 2 112 1 8 /1 fies 
g T7494 61 254 7 mos . 7 3. al 7134 
ce tps 43 260 2 7° 1 686 1.  yalss 
e 7608 62 281 8 108 8 WwW 1 Yeee 
d 6627 58 225 7 106 8 5 2 6501 
Cc 5310 52 190 6 92 — @D 1 ‘6eR2 
me  3f15 35 127 5 85 7 4 | 3508 
a 2139 21 77 4 53 5 2 1 £42011 
5 10,959 105 415 9 103 7 #5 2 10,446 
4 417,003 131 #429704 12 102 7 %5 @ WaeeecO 
3 25,199 159 1316 16 91 7 4A 1 23,795 
2 33,613 183 2438 22 .65 Gf) = ae oT 31,115 
Detector Position 1 < Source Traversed 
Source 
Posit 
11 789 y 33 2 36 4 3 —1 (ies 
10 1346 12 52 2 ? 1261 
9 2079 a 80 dy 1966 
8 3113 32 118 6 2962 
7T 4544 48 175 9 4336 
6 6650 58 261 2 6556 
g 6428 57 2hg 11 6148 
f 6650 58 261 2 6356 
e 6282 56 238 11 6011 
gd 5414 52 212 10 5169 
Cc 4108 nS 161 9 3914 
b 2972 31 TOW, 9 28 32 
a 1670 18 59 3 1578 
5 10,177 1®°(101 389 9 9755 
4h 15,420 124 #670 26 TAR 717 
3 22,895 151 #1287 36 Aig 5 (5 
2 29,172 4170 2421 Ag 26,718 





TABLE VII (Continued) 


BF3 Data = Shutters Installed 


Source Position 11 = Detector Traversed 


Tube 
Posit 
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Op 
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Detector Position 11 - Source Traversed 


Source 


Posit 


I~ 
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30, 907 
28 , 260 


12,740 


1620 
25) 


MW ENO ENONNANANNWO 


D3: 


Zz 
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WEEN FLU £OWwW OF FUIW 
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28,460 
26,868 
18,144 
12, 265 
8064 
ive 
8064 
7612 
6434 
5159 
3508 
2024 
5466 
3729 
ouay 
1521 
860 





TABLE VII (Continued) 


BF3 Data = Shutters Installed 


Source in Position 6f - Detector Traversed 


Tube Np OB No On Np Oo, No OC. Neh Orn 
Posit 
mi} 7610 63 282 5 52 5 2 1 “fis 63 
10 12,197 #=$110 #=«24534 7 78 6 3 1 11,588 #110 
9 23,072 152 Gee 10 86 7 5 2 22,059 4152 
8 35,174 187 1693 13 117 8 4 dW 9333363) Mie. 
7 47,659 218 2933 24 107 8 4 41 £444,623 #4220 
5 48,479 220 2959 24 103 7 5 2 45,422 £«223 
i 36,423 191 1722 #19 #102 7 #=+S 2 Ba scolec 
3 23,346 152 £910 10 91 7 4 1 22,349 #£4x«®152 
2 14,453 120 #496 ve 65 6 4 1 £413,896 120 
1 6548 81 258 5 36 “ 2 | [6257 81 
Detector in Position 6f - Source Traversed 
Source 
Posit 
11 #428112 64 285 12 Whe 1° 5 930) sae 65 
10 14,974 ##122 #4507 22 14,360 124 
9 24,815 157 #908 30 23,800 160 
8 37,077 192 1735 42 35;25p 107 
7 49,947 4223 2940 £54 46,900 229 
5 50,203 224 2897 54 47,199 230 
4 36,122 190 1673 £421 34,342 194 
3 23,024 152 904 30 22,013 154 
2 13,971 #118 #£4518 23 13,346 121 
1 7488 86 258 12 7123 87 


Notes: 
1) Np has been corrected for resolving time. 


(2) Nth is not corrected for diffusion harmonics 
or end leakage 


(3) All data are rounded to the nearest whole count. 
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TABLE VIII 


Counter Stability Tests 


Stability checks were made prior to each complete run, using 


total bare count as the basis for comparison. 


and unloaded pile, without shutters, 
with the source in the pedestal, the detector at station 6(f); 
in the loaded pile with shutters installed, stability was 
checked with the source at 6(f), the detector at 1: 


CeO AY av ewn e bE 


es ee ee ee 
wm £F WwW WH F O 


Loaded, 


Condition 
Loaded 


shutters installed 


DD - 


Total counts per minute 


2027 £18 
2043 £14 
2038 #14 
2047 £15 
6548 £81 
6629 £82 
6683 £82 
6533 £81 
6650 £58 
6568 #81 
6521 £81 
851 £9 
874 £12 
881 £12 
858 £12 


In the loaded 
Stability was checked 


bare 





APPENDIX II 
STATISTICAL EVALUATION OF RESULTS 


The purpose of this appendix is to discuss in some 
detail the statistical treatment of data taken during the 
experiment and to justify the conclusions by means of 


recognized statistical methods. 


A. GENERAL 

As previously mentioned, an attempt was made to achieve 
1% counting statistics throughout the series of measurements. 
In certain cases, however, a total count of 10,000 would 
have required a prohibitive amount of time. In addition 
to variation in counting statistics, it was necessary to 
consider the unequal weighting which must be applied to 
equally valid points on a logarithmic plot. It is desired 
to find the relationship between standard error @& and 
weighting factor w such that the squared deviation of the 
weighted mean is minimized. This is the common problem of 
the constrained extreme, and it can easily be shown that 
with 


Sowpen 


iL 


that the squared deviation is minimized in general when 


K 


Sor 





Now since the standard error of a logarithm is just the 
fractional standard error of the antilog (d(log.x) = 9% 


it follows that 


which is identically the inverse of the variance in the 
logarithm of x. 
The standard error, or standard deviation of the mean, 


has been used throughout this thesis in presenting results. 


B. DETERMINATION OF SLOPES 

In the least-squares curve fitting procedure, stations 
3 through 11 were used in calculating the values of para- 
meters summarized in Appendix I. This was done to allow 
a direct comparison of numbers obtained by the two methods. 
Of the four slopes used as bases for computation, however, 
only the data taken on the loaded pile by conventional 
methods showed a marked deviation from linear regression. 
This occurred at stations 1 and 2, where the detector was 
pce ted within the slowing-down region and near a reflecting 
interface. The fact that proximity of source to detector 
Gid not affect the linearity of the other measurements 
was verified by least-squares analyses using all points. 
In these three cases, no significant change in slope 
resulted. 


The equation used to find the regression coefficient 


Dis 





for minimum squared error is 


24 2494 - nzy 
|| (ee 
2W4242 - ne 


where m 


Slope in linear relation log. (Nip,) @ y = mz+c 


= 
i 


weighting factor 
n= Dw, 
52 2425 
n 
re 
n 
The lines shown in figures 6 and 9 were constructed 


by passing a line of slope m through the intersection of 


xX and jy. 
C. ERROR 


The error in slope was found taking the partial deri- 
vative of m with respect to y, using the standard error in 
the logarithm as the error in y. 

The error in regression coefficient found by this 
method clearly is a function of the counting statistics 
alone, and the slightly smaller error quoted for the recipro- 
cal method reflects only the increase in counting rate 
already mentioned. A similar error could be obtained for 
either method with sufficient time. Much more significant 
is the residual, or the square of the difference between an. 
observed point and the corresponding point Breanicved by the 


least-squares slope. 
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The weighted mean squared errors are tabulated here 


for comparison: 


> Wy (10g,.(N,)) observed ~- log.(Ntp) calculated) ,@ 


2 Wy 


Conventional Reciprocal 
Loaded 4.50x1074 1.27x1074 
Unloaded 3.27x1074 1.30x1074 


D. CORRELATION OF DATA 

An additional value which relates to precision of 
measurement is the value of the slope itself. The student's 
"¢" test, while not devised for this exact application, will 
give a measure of the degree of correlation of the points 


with a postulated straight line (G-2). By definition, 


Mm 
Ff Dey); 
(F-2)2(2;-2)* 


Sy = frog. (Nin) observed ~ loge(Ntn) calculated] 





where 


and F is the number of degrees of freedom, 9 in this case. 


Values of "t" obtained from the data are tabulated for com- 





parison. 
Degree of Data Correlation (t) 
Conventional Reciprocal 
Loaded 82.4 201.5 
Unloaded 141.4 212.0 
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One should note that the statistical parameters mentioned 
thus far can be somewhat BiG soy The standard error of 
the slope is a function only of the counting statistics of 
the points, and will not be altered by poor correlation. 

The mean squared error and the "t" test, on the other hand, 
are functions only of correlation and contain no informa- 
tion on the validity of the points. Considerable research 
has disclosed no parameter which is a function of both of 
these variables. Therefore the authors propose the following 
test of correlation: 

1. Assume a normal distribution of measurements of 

the mean exists about a "true" mean value pre- 


dicted for a point by a least-squares regression 
coefficient. 


2. Define (2) =4 La 


where uy; is the ith deviation from the predicted 
value, O is the standard deviation (or error) 
determined by total count, and n is the number of 
points considered. 


3. Evans (E-1,p.750), shows the probability of ex- 
ceeding any given u/o in a normal distribution. 
From this, the data correlation may be judged by 
the probability that (u/G@ ) will be exceeded in 
subsequent experiments. For the loaded pile, the 
dispersion of conventional data shows a probability 
of being exceeded of only about 15%. More disper- 
sion could be expected from reciprocal data in 37% 
of further experiments. A value near 50% would be 
expected of ideal precision. . 


E. SIGNIFICANCE OF DATA 
In justifying the acceptance of reciprocal data from the 
two lower stations (pile loaded), a normal distribution of 


experimental results is assumed about a value predicted by 
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the slope and the intersection of weighted mean values. 


Calculated Value Observed Value 
Sta. 2 9.3080 * .0085 Si SALES. 
Stm. i 9.6159 + .0070 9.6118 


Sta. 2: u/or = .34: Pu* = .73 
indicating that 73% of additional measurements of log.(N¢n) 
at station 2 would fall farther from the calculated mean, 


27% closer than value observed. 


Sta. 1: u/fo- = .585; Pu = .56 


with similar interpretations. 


F, CORRELATION OF RESULTS 
Results of the experiment are compared similarly, 
assuming: 


1. The value obtained by the conventional 
method is the "correct" value 


2. Additional measurements by the same method 


would follow a normal distribution about 
this "correct" value. 


Conventional Reciprocal ufo Pu 
Material Buckling (cm72) 86 + 12x107° 80 +4 11x107-© .50 .62 


Diffusion Length (em), 
Unloaded 438.7 £ 1.0 48.1 £ .9 60D 





*Pu is defined as the probability that in subsequent measure- 
ments of the same value, by means of identical techniques 
a value in excess of the previously measured value should 
be obtained. Thus values of Pu are generally expected to 
fall within the range .90 to .10. The source of most of 
this material, including notation, is (E-1,Chapters 26-28). 


on i 





Specifically, these values of Pu indicate that the 
parameters obtained by the reciprocal method might have been 


obtained simply by repeating the conventional experiment. 


G. SUMMARY 
In conclusion, this appendix has shown statistical 
proof that: 


1. The pile parameters measured by the two methods 
are statistically identical. 


2. All points of reference in the pile provide 
equally valid data except that in the loaded 
case, data points by conventional measurement 
within 100 cm of the source do not indicate 
linear regression. 


3 The reciprocal method, with similar statistical 


accuracy, provides markedly greater experimental 
precision. i 
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APPENDIX III 
DETAILS OF PROCEDURE 


This Appendix discusses in detail the experimental 


methods and calculation procedures used for indium foil 


measurements. 
1. Equipment 


Indium foils of 7/8" diameter were punched from .005" 
stock. Cadmium covers used in measuring the epithermal 
flux component were 1" in diameter and .020" thick. In 
order to minimize effects of counter background, beta radia- 
tion from the foils was employed as a measure of activation. 
Counting was done with a proportional flow counter of 27 
geometry. The counter is manufactured by Nuclear Chicago and 
uses a 90% argon, 10% methane gas. This counting method 
possesses advantages of high efficiency, stable geometry, and 
excellent time resolution. The counter was connected to 
an automatic sample changing mechanism which was preset for 
1000 or 400 counts, depending upon the expected degree of 
foil activation. Counter stability was checked frequently 
by measurement of a standard laboratory beta source. Suffi- 
cient counter background measurements were made to establish 


this factor with essentially no error. 


2. Calculation Methods 
A printed tape received from the changing and recording 


mechanism provides the exact time required for a preset number 


ote 








of counts from each foil. Elapsed time since removal from 
flux is easily obtained from this, knowing precisely the 
time required for sample changing. If x' is the time ans 
quired to obtain a given number of counts with background 
present and z is the time required for the same number of 
counts from counter background alone, it can easily be 


shown that 


where x is the time which would be required for this number 
of counts without background. 

Since the neutron flux desired is proportional to 
saturation foil activity, it is necessary to correct further 
to allow for the following: | 


a, Lack of complete saturation activity 
due to inadequate exposure time. 


b. Decay of radioactive isotope during 
the period following removal from flux 
and prior to counting. 
c. Decay during the counting interval. 
The mechanics of these corrections are well known and do 
not require discussion. For measurements reported in this 
thesis, the 13 second isotope of indium was allowed to decay 
before counting and only the 54.2 minute activity was measured. 


A form of the decay equation suitable for repetitive 


calculation is: 
N Xo PEL e Aex 
xX 
pene 





Ao da = 
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wheres 
AjAq = saturation activity (CPM) 


N = preset counts 

Aa = decay constant of 54.2 minute 
isotope 

ty zs @lapsed time from foil removal 


to start of counting period 


x = counting time without background 


This expression assumes that exposure time was sufficient to 
obtain saturation of the isotope being measured. Such was 
the case for all data included in this report. 

After corrected saturation activities were obtained, 
Standard methods were used to reduce these data to relative 
neutron flux density. A convenient method for this is 


reported by Wilson, et al (W-3): 


where cadmium ratio is measured with source in position. The 
above expression simply assumes that CR does not change in the 
no-source condition. This supposition is not valid, but its 
use avoids cadmium covered measurements of the spontaneous 
fission pile background. These generally are not very satis- 
fying from the statistical point of view. The simplification 
is justified by the fact that error introduced thereby is very 


small. 


3. Results 


Counting rates from indium decay induced by thermal neutrons 
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varied from 12.5% to 71.7% of the total registered count. Hence 
statistical variation in the measurements rendered them useless 
for accurate calculation of "control" parameters. It was not 
considered worthwhile to correct foil results for slowing down 
effects, diffusion harmonics, or end leakage. The data, as 
reported in Table IV, are corrected to relative thermal neutron 
density. Foils were counted front and back to eliminate 
shadowing effects. Certain refinements (D-2,M-3,T-1) which 
might have been applied were not statistically justified. No 
attempt was made to calculate pile parameters from foil measure- 


ments. 
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APPENDIX IV 
SUPPLEMENTARY DISCUSSION 


The purpose of this appendix is to discuss qualitatively 
Some aspects of the experiment which do not bear directly upon 
the objectives of this thesis. ~_ 

several features of the experiments conducted by cone 
ventional and novel means warrant consideration. There is a 
Separation of lines in figure 9, where thermal neutron density 
in the unloaded pile measured by the reciprocal method was 
about 10% greater than that measured by the conventional 
method. Second, there is a separation of lines in figure 6, 
showing that in the loaded pile, an increase of about 77% of 
conventional thermal neutron count rate was observed in the 
reciprocal experiment. In addition, the conventional method 
produced a disproportionate increase in counting rate at 
stations 1 and 2 in the loaded pile which was not observed 
using the reciprocal method. Third, there is a depression 
at the center of the conventional horizontal traverses (figures 
5 and 8) which was not observed by the reciprocal method. 
Finally, with both source and detector in the lattice and with 
cadmium shutters installed, there is a difference in count 
rate between the two methods, which reverses as the reference 
point is changed from station 1 to station 11. The authors 
have found no rigorous explanation of the latter two phenomena. 
In relation to the deviation of horizontal distributions 


from the expected cosine, the authors have (1) investigated 


67. 





the pile for inclusions, slots, irregularities of any des- 
cription (2) attempted fitting curve with harmonics of the 
fundamental mode of distribution (3) rechecked raw data. None 
of these steps have revealed the cause of the deviation. | 

The first of these conditions is ascribed to differences 
in net neutron loss from the system. Such an effect would 
be expected where the system is not isolated from its 
Surroundings. The pedestal rests on a felt base on the 
concrete foundation; there is no cadmium between the graphite 
and the foundation of the building. Taking this 10% difference 
in neutron density as a function of geometry and surroundings, 
we would expect the factor to apply similarly to the loaded 
assembly. 

As shown by figure 8, the separation of actual count 
rates is much greater in the loaded assembly than the geo- 
metry factor would predict, i.e., the apparent multiplication 
of neutrons by loading natural uranium in the lattice is 
much greater in the reciprocal case than in the conventional 
case(see Appendix V). This is not predicted by one-group 
theory, which assumes a thermal source. Time limitations 
have prevented the authors from performing multigroup calcu- 
lations which might be expected to disclose the reasons for 
the increased count. 

Part of the increase may be attributed to fast fission of 
u238 caused by source neutrons. This could not be expected to 


exceed 5% according to calculations which conservatively 
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aSsume a fast flux equal to thermal flux. 

The streaming effect is not considered to produce a em 
nificant difference between loaded conventional and reciprocal 
counts. The pedestal channel extends through the pile, while 
the detector channels extend only about 3/4 of the horizontal 
dimension. The Simon-Clifford equation (GR-3,p.276) gives 
a ratio of streaming to total flux of only about 5 x 107, 
indicating a negligible difference in streaming whether the 
source is in the pedestal or the lattice. 

The great difference in cadmium ratio between conven- 
tionally and reciprocally measured neutron densities, shown 
in figure 10, tells us that the average neutron energy is high 
in the lattice, diminishes near the interface, and is near 
thermal in the pedestal. This leads to investigation of 
average neutron temperature in the two regions. Assuming 
an equal total neutron population with Naxwell-Bolt2natin 
distribution in both cases, and assuming, from unloaded data, 
al.l geometrical advantage in the reciprocal case, the re- 
ciprocal total count is 50% higher than expected. Suppose this 
is a result of difference in average energy seen by the 1/v 
boron. If the higher count reflected a kT of .025 ev, the 
lower value would mean a cross section of about .67 of 
"thermal" value, hence a kT of about .056 ev. This corresponds 
to a neutron temperature of about 380°C, which is not unusual 
for a reproducing region. It is a mistake, however, to attach 


too much significance to numbers produced by the assumption 
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STATIONS ON Z-AXIS 





of Maxwell-Boltzmann distribution. Many experimenters have 
attempted to define the very significant deviation of modera- 
ted neutrons from this distribution (W-1, p.332, and W-2). There 
is poor agreement between theory and experiment, and even among 
experiments. | 

A difference in flux "hardness" in the vicinity of the 
detector, then, is supported qualitatively by theory. Phy- 
Ssically, we should expect the average neutron energy to be 
greater in a reproducing region than in moderator. Further- 
more, the flux would tend to decrease in average energy near 
a reflecting interface resulting in the trends observed in 
figures 6 and 10. The authors are of the opinion that the 
shift in average neutron energy, evidenced by cadmium ratio, 
is the primary reason for the differences in counting rates 
by the two methods. It is cautiously suggested that experi- 
ments of this type might be designed to study neutron spectra 


in lattices. 
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APPENDIX V 
MULTIPLICATION RATIO 


The exponential experiment is only one of two basic 
methods for determination of the critical size of a reactor. 
The other procedure, which will be briefly considered here, 
is termed the "critical assembly" method. This method requires 
the construction of a flexible fuel-moderator assembly which 
can be gradually increased in size until a self-sustaining 
reaction is attained. Criticality is predicted by plotting 
the inverse multiplication ratio versus fuel mass as the load- 
ing is successively increased. The fuel mass corresponding 
to a zero value of reciprocal multiplication ratio (infinite 
multiplication) indicates the fuel requirement for criticality. 

such a procedure is not feasible with a fixed subcritical 
assembly of the type used during this thesis. However, it is 
felt that values of multiplication ratio for the single fixed 
configuration are of interest and should be reported. 

This multiplication can be determined experimentally by 
measurement of the thermal neutron density within an assembly 
at a specified distance from the external source. This measure- 
ment is made at each location with the pile loaded and sub- 
sequently after all fuel has been removed. The ratio of 
these measurements gives the required multiplication, ata 
fixed point, for the particular size of the assembly which 
was investigated. 

Such measurements were made by conventional and reciprocal 


ee 





a 2 


methods for the M.I.T. subcritical graphite pile. Values 
were obtained for this ratio of loaded to unloaded thermal 
neutron density at each of the measurement positions employed 
throughout the thesis. Figure 11 illustrates the results 


which were obtained. 
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APPENDIX VI 
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summary of information published prior to 1959 on the sub- 
ject of exponential experiments. Certain papers are included 
which relate to reciprocity and the importance function. By 
no means are all of these works cited within the text of this 
thesis. It is felt, however, that a thorough compilation may 


be of value to others contemplating work in these fields. 
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